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The determination of precision and accuracy in horizontal geodetic networks using 
a least squares adjustment is a complex process that has a range of available 
software applications. However, it is commonly observed that users may not 
possess a complete understanding of the underlying theories and equations that 
inform the development of such software. This paper presents a set of step-by-step 
procedures for determining the precision and accuracy of horizontal geodetic 
networks using the observation equation method of the least squares technique. In 
the sample calculation, the a posteriori variance and standard error were computed 
to obtain the precision and accuracy of the network, respectively. The inverse of the 
normal matrix was multiplied with the a posteriori variance to obtain the precision of 
the adjusted parameter. The square root of the precision was evaluated to obtain 
the accuracy of the adjusted parameter. The enumerated procedures are presented 
in a sequential order to ensure ease of comprehension for users. Additionally, a 
numerical example is provided to illustrate the procedures and to enhance 
understanding of the concepts presented. The computed adjusted network's 
precision and accuracy are 0.723m? and 0.850m respectively. Our objective is to 
provide users with the necessary tools and knowledge required to confidently 
perform this task. 
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The Least Squares technique is an essential tool 
for analyzing and correcting random errors in 
survey measurements. According to Zulkifli et al. 
[1], the Least Squares Adjustment (LSA) method 
is an advanced adjustment technique that refines 
observations based on _ the principles of 
probability, which account for the occurrence of 
random errors. By utilizing estimated precisions 
of observation coordinates, LSA _ calculates 
adjusted positions that reconcile discrepancies 
between the observed and adjusted coordinates. 


Furthermore, this method produces 
comprehensive _ statistical data on _ the 
adjustments made, enabling an assessment of 
the strength and confidence in computed 
positions, while also facilitating the detection of 
potential blunders. As a rigorous approach, it is 
particularly valuable in the adjustment of 
horizontal geodetic networks, yielding the most 
accurate and reliable values for survey 
measurements. The Least Squares technique 
guarantees the determination of the most 
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probable values of observations based on a 
given set of redundant measurements [2]. 


A horizontal geodetic network consists of a 
series of reference stations, each of which has its 
position and coordinates established with a high 
degree of accuracy. The network is categorized 
into three groups: first-order, second-order, and 
third-order. The classification is based on the 
precision with which the networks. are 
established. The accuracy and precision of the 
network are typically calculated using the least 
squares adjustment method. 


Researchers have encountered challenges in 
comprehending the accuracy and precision of the 
horizontal geodetic network when employing the 
least squares technique. However, fortunately, 
there are software programs available to carry 
out the adjustment. Despite this, many users still 
struggle to grasp the theory and equations 
employed in developing these programs. The 
matrix nature of this technique has made it 
particularly challenging. The least squares 
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adjustment can be represented as the sum of the 
estimate matrix and the observation matrix, 
equating to the residual matrix. Unfortunately, 
obtaining each of these matrices has posed a 
significant problem. Previous studies [3,4,5,6] 
that have employed the least squares adjustment 
technique have failed to provide a breakdown of 
how this technique is applied. To assess the 
accuracy and precision of the adjusted network, 
it is necessary to obtain the residual matrix and 
have redundant observations, as well as a 
degree of freedom. 


This paper describes the observation equation 
method of the least squares adjustment 
technique for determining the precision and 
accuracy of horizontal geodetic networks 
established using DGPS. The application is 
presented in a detailed, step-by-step manner. 


1.1. Observation Equation Technique of 
Least Squares Adjustment 


The observation equations are mathematical 
expressions that link observed quantities to 
independent unknown parameters and 
observational residuals. Each observation has a 
unique set of unknowns and requires a 
corresponding equation. To arrive at a unique 
solution for these unknowns, the number of 
equations must equal the number of unknowns. 
In most cases, there are more observations than 
unknowns, allowing the principle of least squares 
to determine the most probable values for the 
unknowns [7]. 


Ayeni [8] and Okwuashi and Asuquo [9] have 
described the observation equation method, 
which expresses adjusted observations as a 
function of adjusted parameters. As outlined in 
[10], the functional relationship between adjusted 
observations and adjusted parameters is 
established through this approach (Equation, 1). 


L,=F(X,) (1) 


Where, 
L,, = adjusted observations 


X , = adjusted parameters 
Equation (1) is a linear function and we have 
derived the general observation equation model. 


To create a matrix expression for conducting a 
least squares adjustment, comparisons will be 


CaJoST, 2024, 3, 265-272 


S. O. Eteje et al. 


draw with systematic procedures. The matrix 
notation for the observation equation system can 
be found in [11] and [12]. 


V=AX-L (2) 


Where, 
A = Design Matrix 
X = Vector of Unknowns 


L = Calculated Values (/,) Minus Observed 
Values (1, ) 
V = Residual Matrix. 


That is, 
es a, Ay Qi, 
V= oa we: Ay, Ano Aon |5 
Vin an Ano Ann 
Xx 
1 L, 
Xx 
X=| °° and yp l, 
Xin Ls 
The process of determining the unknown 


parameters matrix, X involves the deduction of 
the normal matrix, N and the matrix of numeric 
terms, t When dealing with weighted 
observations, the weight matrix, which is always 
a square matrix, also needs to be derived. 


As stated by Ghilani [13], a system of weighted 
linear observation equations can be presented in 
matrix form as: 


A™WAX = ATWL (3) 


To make X the subject of the formula, both sides 
of equation (3) will be divided by A’WA. Thus, 


X =(A'WA)'A'WL (4) 


If A7WA = N, normal matrix and A’WL =1, matrix 
of numeric terms, then equation (4) becomes 


XL =N 4G (5) 
1.2 Weight of Uncorrelated Observation 


Deakin [14] posits that the weight, w of an 
uncorrelated observation is directly proportional 


. . 2 
to the inverse of variance, o~ . Hence, 
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yee (6) 


The weight matrix of uncorrelated observations is 
a diagonal matrix such that the off-diagonal 
elements are all zero and it is presented as 


— 0 0 0 
at 
_| 0 ae 0 O 
O17 
0 0 0 
0 0 O ev 
= 


1.3. Precision and Accuracy Determination 


When calculating the mean precision and 
accuracy of a set of adjusted observations, a 
posteriori variance and standard error are used. 
It's worth noting that variance and standard error 
are measures of precision and accuracy, 
respectively. The formulae for calculating a 
posteriori variance, as outlined in [15] (Equation, 
8), are utilized in this process. 


a. _V'WV 


oO 


(8) 


r 


The model for the a posteriori standard error 
given in [16] and [12] is 


T 
ep aks (9) 


r 


Where, 

r=n-m-= degrees of freedom 

n= number of observations 

m= number of unknown parameters. 


The precision and accuracy of the adjusted 
parameter and coordinates are derived from the 
variance-covariance matrix. The precision of the 
adjusted parameters is represented by the 
variances and elements of the principal diagonal, 
whereas the off-diagonal elements indicate the 
co-variances between specific coordinates. The 
square root of the variances determines the 
accuracy of the adjusted parameters. To 
calculate the variance-covariance matrix, the 
inverse of the normal matrix is multiplied by the a 
posteriori variance. 


= ON (10) 
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The trace is a commonly used metric for 
assessing accuracy in_ statistical models. 
Specifically, it refers to the sum of the diagonal 
elements of the variance-covariance matrix. 


pee It is usually written as rp, or 


sometimes tr(}’.. ), that is, 


CO = 56,0. 


As stated in [17], one can interpret the trace of 
the variance-covariance matrix as a means of 
determining the accuracy of the related set of 
random variables. The calculations for 
determining the variance and standard error of 
the adjusted coordinates can be found in the 
models presented in [12]. 


(11) 


o 1 = Oo, On (1 2) 
Og = ¥F5Onn (13) 
Where, 


Q,,, = inverse of the normal matrix main diagonal 


element. 


Error ellipses are a common calculation in the 
adjustment of a horizontal or three-dimensional 
network. They provide a helpful measurement for 
interpreting the accuracy of positional station 
directions. Error ellipses are also commonly 
applied in optimizing a given network. 

Mikhail and Gracie [18] gave the semi-major axis 


2 oe ee) ; 
1, semi-minor axis o, and the orientation of 


the error ellipse 0 as 


2 4 » (14) 
2 2 2 2\2 
gi, POs, || LOL Her) é 
ar : | Fi —|+o,, 
20 
C,- 0, 
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2. Methodology 


2.1 Steps to Compute Precision and 


Accuracy 


The following steps are necessary for the 
adjustment of a GPS network: 


1. At least two fixed controls must be used to 
observe each network point. The data 
acquired with respect to any of the 
controls/base stations should be processed 
separately to obtain the coordinates of the 
new points. 


2. After obtaining the coordinates of the new 
points, the observation equations must be 


deduced. The number of observation 
equations must not be less than the number 
of baselines. 


3. From the deduced observation equations, 
the coefficient and design matrix (A), 
observation matrix (L), residual matrix (v), 
and matrix of unknown parameters (X) must 
be derived. 


4. The weight matrix should also be deduced 
using Equation (6). The variances of the 
baseline vectors for DGNSS observation are 
obtained from the variance-covariance matrix 
of the processed DGNSS data. 


5. After deducing the matrices, the unknown 
parameters must be computed using 
Equation (5). 


6. The design matrix is used to multiply the 
computed values of the unknown parameters 
to determine the estimate as well as the 
most probable values. 


7. The residuals must be evaluated by finding 
the differences between the estimates and 
the observations using Equation (2). 
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8. The precision of the adjusted observations 
can be determined using Equation (8) after 
computing the residuals. The degree of 
freedom must be calculated to obtain the 
precision of the adjusted observations. 


9. The mean accuracy of the adjusted 
observations can be determined using 
Equation (9) by taking the square root of the 
determined precision. 


10. The accuracy of the adjusted network can 


also be determined using Equation (11). 


11. The precision and accuracy of the adjusted 
coordinates can be determined using 
Equations (12) and (13), respectively, after 
determining the precision and accuracy of 
the adjusted observations. 


12. The semi-major axis, semi-minor axis, and 
orientation of the error ellipse can be 
computed using Equations (14) and (15), 


respectively. 
3. Results and Discussion 


3.1. A Sample Calculation of Precision and 


Accuracy 


Table 1 shows the rectangular coordinates and 
corresponding variances for network stations A 
and B in relation to controls S and T, as 
determined through the use of DGPS. The 
positions for controls S and T_ remain 
350472.960mE, 251374.548mN, and 
354095.611mE, 251441.978mN correspondingly. 
Compute the most probable positions for points 
A and B, as well as the accuracy and precision of 
the observations and those of the adjusted 
positions using the least squares method. 
Additionally, compute the error ellipses for the 
coordinates of the stations. 


Table 1: Observed Positions and their Respective Variances 


spate | FoNer] Coordinatesim) | Vaianee | Base | Rover | coordinates im) | Varnes 
is Northing | 950852.942 | 0.0000577 , Northing | 250852.949 | 0.0000465 

é Easting | 352598.178 | 0.0000314 7 Easting | 352598.188 | 0.0000554 
Northing | 252127.392 | 0.0000247 ; Northing | 252127.398 | 0.0000338 

Easting | 352572.216 | 0.0000822 Easting | 352572.226 | 0.0000709 
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Figure 1: Observed Network Points 


three-dimensional coordinates from a redundant 
set of measurements. 


The method of least squares adjustment readily 
satisfies the requirement to accurately estimate 


Solution 
Table 2: Determination of baseline vector 
COORDINATES 
STATION STATION 
FROM AN (m) AE (m) NORTHING (m) EASTING (m) TO 
251374.548 350472.960 s 
Ss -521.606 2125.218 250852.942 352598. 178 A 
Ss 752.844 2099.256 252127.392 352572.216 B 
251441.978 354095.611 T 
T -589.029 -1497.423 250852.949 352598. 188 A 
T 685.420 -1523.385 252127.398 352572.226 B 
Deduction of observation equations 
Ay =Sy +AN,, +V, 17331.02 0 0 0 0 0 0 
0 3184713 0 0 0 0 0 0 
gs ale a 0 0 — 40485.83 0 0 0 0 
By =Sy + ANsy +V; wi ° 0 0 12165450 0 0 
B, =S, +AEg, +V, “| 0 0 0 0 2150538 =O 0 
0 0 0 0 0 18050.54 0 
Ay =Ty + AN yg + V5 0 0 0 0 0 0  —-29585.80 0. 
A, =T, + AE,, +V¢ 0 0 0 0 0 0 14104.37 
By =Ty +ANi, + V; Ay 
B, =T, + AE, + Vz : A; 
|B 
Formulation of observation matrix, L, design s 
matrix, A, residual matrix, v, weight matrix, w and e 
matrix of unknown parameters, X, ; 
7 Computation of the unknown parameters, X 
Ay 0 0 0 1000 using Equation (5) 
0 A, O O 0 1 00 38836.40 0 0 0 
By 8 ce’ a Normal matrix, N = A’WA = : ae ? i" 
0 0 0 B, 0001 0 0 70071.63 0 
A= = 
A 0 0 0 1000 0 0 0 26269.82 
0 A, O O 0 1 0 0 , 
0 O By Of} }0 0 1 0 0.0000257 0 0 0 
0 O O B, 0001 3 ee 0 0.0000200 0 0 
N =(A'WA)" = 
Sy + AN) (250852.942 Yi : i ae 
x “ : Vv. 0 0 0 0.000038 1 
S, + AEs, 352598.178 2 
V. 
BON ae | | Ae ee oh 9742225063.763 
S, t+ AEs, | | 352572.216 y.| "4 
pe = ella 7 17593829570.669 
Ty + AN, 250852.949 5 Matrix of numeric terms, t = A’ WL = 
7, + Ag, | |352508.188 V, 17666977195.008 
T, +AN,, | | 252127.398 V, 9262009666.744 
T, + AE np 352572.226 V, 
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0.0000257 0 0 9742225063.763 250852.946 
~yeNn+4ye 0 0.0000200 0 175938295 70.669 2 352598.182 
0 0 0.0000143 17666977195.008 252127.395 
0 0 0 0.0000381 9262009666.744 352572.221 
0.0000257 0 0 0 
Calculation of the most probable positions for 0 0.0000200 0 0 
stations A and B Dive = 0.723 0 0 0.0000143 0 
1 00 0 250852.946 0 0 0 0.0000381 
0 1 0 0 352598.182 
0 0 1 0 250852.946 252127.395 0.0000186 0 0 0 
AX = 00 0 1 ' 352598.182 a 352572.221 2 0 0.0000145 0 0 
100 0 252127.395 250852.946 0 0 0.0000103 0 
0 1 0 O| (352572.221} | 352598.182 o : 0 GONE 
De Get E eal Thus, the accuracy of the adjusted network using 
00 0 1 352572.221 


Therefore, 250852.946mN, 352598.182mE and 
252127.395mN, 352572.221mE are the adjusted 
as well as the most probable positions for points 
A and B respectively. 


Evaluation of residuals, v using Equation (2) 


250852.946 250852.942 0.004 
352598.182 352598.178 0.004 
252127.395 252127.392 0.003 
yay p= 352572.221 | | 352572.216 i: 0.005 
250852.946 250852.949 —0.003 
352598.182 352598.188 — 0.006 
252127.395 252127.398 —0.003 
352572.221 352572.226 — 0.0005 


The precision and accuracy of the adjusted 
network can be computed using Equations (8) 
and (9) respectively. As there are 4 unknowns 
and 8 observations, the degree of freedom, 
denoted by r, is equal to 4. 


a= 2.891 = 0.723m? 
4 
- = = 40.723 = 0.850m 


The adjusted network's precision and accuracy 
are 0.723m? and 0.850m, respectively. The 
accuracy of the adjusted network was computed 
using Equation (11) which is the trace of the 
variance-covariance matrix. 
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Equation (10) is 
y = 0.0000186 + 0.0000145 + 0.0000103 
XX +.0.0000275 = 0.0000709m7 


=~/0.0000709 =0.00842m 


Calculation of the precision and accuracy of the 
adjusted coordinates using Equations (12) and 
(13). 


. = Ay 0.0000257 0.0000186 
Station A precision, = 0.723} = 


Ap 0.0000200) |.0.0000145 
; = " erie aes 
Station B precision, = 0.723 = 
B, 0.0000381) \.0.0000275 
: Ay 0.0000186 0.00431 
Station A accuracy, = = 
A, ) \0.0000145) (0.00381 
Station B accuracy, 
& _ ({0.0000103 | (0.00321 
B, ) \0.0000275 ) \0.00524 


Hence, the precision of northing and easting 
coordinates of stations A and B are respectively 
0.0000186m2, 0.0000145m? and 0.0000103m?, 
0.0000275m? while their respective accuracy are 
0.00431m, 0.00381m and 0.00321m, 0.00524m. 


Calculation of error 
positions using 
respectively. 


ellipses of the station 
equations (14) and (15) 


For station A 
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>  0.0000186+0.0000145 | (0.0000186 —0.0000145)" 
o1= - + ri +0 


= 0.000033m? 
es 0.0000186 + 0.0000145 
y 2 
= 0.000033m? 
o , = ¥0.000033 = 0.00575m 


o@ , = V0.000033 = 0.00575m 


0 
>6=0° 
0.0000186 —0.0000145 


| coonore 0000s ais 
4 


tan 20 = 


The procedure was carried out again for Station 
B and the following values were obtained. 


o , =0.00615m 
o, = 0.00615m 
0=0° 


Accordingly, the semi-major axis, semi-minor 
axis and the orientation of positions A and B are 


0.00575m, 0.00575m, 0° 
0.00615m, 0.00615m, 0°, respectively. 


and 


It should be noted that there is no covariance 
between stations A and B due to _ the 
uncorrelated nature of the observations. This 
results in zero orientation of the error ellipses. 


4. Conclusion 


The precision and accuracy of a horizontal geodetic 
network can be determined through the rigorous least 
squares method, which also provides information on 
the network's reliability and classification. This paper 
offers a step-by-step guide to applying the rigorous 
method and showcases a numeric example for a 
better understanding of the processing procedures. By 
following these procedures and examining the 
provided example, users can gain a comprehensive 
understanding of the least squares adjustment 
software and its associated theories and equations. 
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